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Abstract 


One of the objectives of hydrological studies is the estimation of the 
annual runoff volume (VEA, for its name in Spanish) of watersheds 


without gauging. Such data is basicin the design of reservoirs. In this 
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study, the estimation is provided through regional frequency analysis, 
based on the Runoff Index method. This approach was applied in the 
Hydrological Region No. 10 (Sinaloa), Mexico, for which 22 VEA records 
were processed, with amplitudes that varied from 24 to 57 years. 
Through the Discordance Test, three records were eliminated. Regional 
hydrological homogeneity was verified with the remaining 19 records 
and three statistical techniques: seasonality indices, linearity of 
ordinary moments, and linear regression between probability-weighted 
moments. The tests of loss of homogeneity of each record led to the 
elimination of four, due to their ascending linear trend. The 15 
remaining records were processed in a dimensionless way, by dividing 
theirdata over their annual mean runoff volume (VEMA) and they were 
concatenated with the stations-years method, forming a series of 539 
data. Three probabilistic models were applied to obtain the regional 
growth curve: Pearson type III, Log-Normal, and Potential 
Transformation; choosing the one with the smallest standard error of 
fit. The scaling of the regional Log-Normal distribution was achieved 
with a logarithmic relationship between the VEMA and the basin area 
of the processed records. In the Conclusions section, the established 
procedure is detailed and the importance of the results to generate 
synthetic sequences of VEA is highlighted (which are equally probable 
to occur). Therefore, the systematic application of the Runoff Index 
procedure as a method of estimating the VEA in other regions of the 


country is recommended. 
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Resumen 


Un objetivo de los estudios hidrológicos es la estimación del volumen 
escurrido anual (VEA) de cuencas sin aforos. Datos que son básicos en 
el diseño de embalses. En este estudio, tal problema se resuelve por 
medio del análisis de frecuencias regional, con base en el método del 
índice de escurrimientos. Este enfoque se aplicó en la Región 
Hidrológica No. 10 (Sinaloa), México, procesando 22 registros de VEA, 
con amplitudes que variaron de 24 a 57 años. Por medio de la prueba 
de discordancia se eliminaron tres registros. La homogeneidad 
hidrológica regional se verificó con los 19 registros restantes y tres 
técnicas estadísticas: índices de estacionalidad, linealidad de los 
momentos ordinarios y regresión lineal entre momentos de 
probabilidad ponderada. Las pruebas de pérdida de homogeneidad de 
cada registro condujeron a eliminar cuatro, por mostrar tendencia 
lineal ascendente. Los 15 registros remanentes se procesaron de forma 
adimensional, al dividir sus datos entre su volumen escurrido medio 
anual (VEMA) y se concatenaron con el método de las estaciones-años, 
formando una senie de 539 datos. Para obtener la curva de crecimiento 
regional se aplicaron tres modelos probabilísticos: Pearson tipo III, 
Log-Normal y Transformación Potencial; adoptando el de menor error 


estándar de ajuste. El escalamiento de la distribución Log-Normal 
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regional se logró con una relación logarítmica entre el VEMA y el área 
de cuenca de los registros procesados. Las Conclusiones detallan el 
procedimiento establecido y destacan la importancia de los resultados 
para generar secuencias sintéticas de VEA, igualmente probables de 
ocurrir. Por ello, se recomienda la aplicación sistemática del método 
del Índice de escurrimiento en otras regiones del país, para contar con 
tal método de estimación del VEA. 


Palabras clave: volumen escurrido anual, momentos L, prueba de 
discordancia, índices de estacionalidad, homogeneidad hidrológica 
regional, tendencia lineal, error estándar de ajuste, secuencias 
sintéticas. 
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Introduction 


Reservoir sizing 
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Surface Hydrology, about the estimation of hydraulic resources, seeks 
to answer the following two questions: how much water is available? 
and what are the extreme values? Both issues involve estimating the 
available volumes of water and ¡its extreme magnitudes (floods and 


droughts), as well as their occurrence or behavior over time. 


The classic hydraulic works for utilization are the reservoirs, 
which can be divided into large and small; the former is multipurpose 
and generally store water from the wet years and use it in the dry 
years, the latter store water from the rainy seasons and release it in 
the next dry season. The hydrological design of reservoirs, whether 
large or small, requires the estimation ofthe annual and monthly runoff 
regime at the project site, for their sizing. Furthermore, designing 
floods are necessary to guarantee their hydrological safety (Carr € 
Underhill, 1974; Zsuffa € Gálai, 1987; Skertchly-Molina, 1989). 


Due to the low density of hydrometric stations, the runoff regime 
estimation atthe reservoir project site was generally carried out before 
the regional approach, based on the concept of similarity of basins. 
Then, a hydrological model of the rain-runoff process, calibrated in a 
gauged basin, was applied to another similarbut without hydrometric 
data (Cannarozzo, Noto, Viola, 8 La-Loggia, 2009). 


Currently, the regional frequencies analysis allows the use of all 
the hydrometric information available in a geographic area, through 
statistical methods or techniques which help define a homogeneous 


region or zone from a hydrological point of view, and then it is feasible 
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and reliable to establish a common  dimensionless probability 
distribution function (PDF), about a mean or median índex value. 
Subsequently, in each homogeneous region, an empirical relationship 
is obtained between the index value and various morphological and 
climatic characteristics ofthe basin easily calculable, making it possible 
to generate syntheticsamples of the annual runoff volume (VEA, by its 
acronym in Spanish) at any site in the studied region (Cannarozzo et 
al., 2009). Blóschl, Sivapalan, Wagener, Viglione, and Savenije (2013) 
indicate that the hypothesis of this method establishes that the annual 
mean runoff volume can vary between sites in a homogeneous region, 
within which the PDF is identical and defined by the dimensionless 


growth curve. 


On the other hand, Skertchly-Molina (1989) has presented a 
thorough classification of the methods that allow dimensioning the 
reservoirs. Such techniques can be divided into two types: (1) those 
which process historical runoff records and (2) those which use 
synthetic records that are equally likely to occur. Regarding the 
convenience of one or the other, McMahon and Mein (1986) point out 
that the historical record of runoff is the least probable or feasible 
method that will repeat exactly in the future. In addition, they indicate 
that having a way to view the sensitivity of the reservoirís necessary 
capacity estimate is very useful to make more reliable decisions. This 
is possible by making use of syntheticseries of annual runoff (Campos - 
Aranda, 2011). 
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Objective 


The present study focuses on the regional frequency analysis of the 
annual runoff volume (VEA, by its acronym in Spanish), in millions of 
cubic meters (Mm3), from 22 hydrometric stations of the Hydrological 
Region No. 10 (Sinaloa), Mexico. The results allow the generation of 
synthetic sequences, of any size, equally likely to occur from the VEA, 
at the sites of the hydrometric stations used, and at any location 
without gauging data within the hydrological region analyzed. The 
Runoff Index method was used, similarto the well-known Flood Index. 
We worked with 19 records, after applying the Disagreement Test and 
the regional hydrological homogeneity was verified based on three 
statistical techniques. Next, the regional PDF was defined, through the 
standard error of fit, and finally, a logarithmic equation that relates the 


annual mean runoff volume and the basin area was obtained. 


Methods and materials 
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Runoff Index Method 


Hosking and Wallis (1997) mention that the techniques integrating the 
so-called index flood method can be applied to any type of data. This 
method was proposed by Dalrymple (1960) to process flood data, also 
called maximum flows (hence its name) allows estimating a summary 
of statistical parameters from the grouping of several different samples 
or records (Rahman, Haddad, € Eslamian, 2014; Ouarda, 2017). 
When there are records of annual runoff volumes at NS sites and 
each location j, the number of data is ni of volumes ve, Then, we have 
that V¡(F) with O<F<1, is the quantile function of the PDF of each site 
j. The basic hypothesis of the runoff index method establishes that the 
sites j integrate a homogeneous region, which implies that the PDFs of 


the NS sites are identical, except for the scaling factor, which is 
precisely the Runoff Index (IE). The equation of the method is: 


VI(E) = IE! -q(F) with j = 1,2,...,NS (1) 


where q(fF) is the regional growth curve, that is, a dimensionless 
quantile function common to all sites, and q(F) is a regional PDF or 


common function of the values V?/1E), In general, the IP is equal to 
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the mean annual runoff volume (VEMA/, by its acronym in Spanish) of 
the data from each site, but the median can be used (Hosking € Wallis, 
1997). 


In the 7£/ method, the following four hypotheses are accepted: 
(1) the runoff from each site are independent, (2) the runoff from 
different sites is independent, and (3) the PDFs at different sites are 
identical except for the scaling factor and (4) the mathematical form 
of the regional growth curve is perfectly specified (Hosking 8 Wallis, 
1997). 


The regional frequencies analysis, including the IE, has two 
fundamental applications, within the homogeneous region; the first is 
carried out in sites with short or unreliable records and the second in 
localities without hydrometric data, using the empirical equations 
obtained for the relationship between the VEMA/ and the basin area, or 
other morphological characteristics of the basin, or the main channel. 


Moments and L ratios 


The L moments (A;) are linear combinations of the weighted probability 
moments (fr) developed by Greenwood, Landwehr, Matalas, and Wallis 


(1979), which are statistical parameters of the ordered data. The L 
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moments have proven to be an efficient and robust system for the fit 
of PDFs currently used in hydrology or established under the precept. 
They are calculated with the following equations (Hosking 8 Wallis, 
1997: Rao € Hamed, 2000; Stedinger, 2017): 


M= Bo (2) 
%. = 2" B.—Po (3) 
A¿=6-P2—6-B,+Bo (4) 
Ay =20-B¿— 30-B, +12: B,—Bo (5) 


In such a system, the quotients (7) of L moments are defined, 
starting with L-Cv which is analogous tothis coefficient, and then those 
of similarity with the asymmetry and kurtosis coefficients, these are: 


To =2%/A; (6) 
T¿ =23/A, (7) 
T, =2y/A, (8) 
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In a sample of size n, with its elements X; exposed in ascending 
order (X, < X, < *** <X, ) the unbiased estimators of B-are obtained with 


the following general expression: 


ly E _ 
B, = A a a X, conr = 0,1,...4 (9) 


Discordancy test 


In any regional frequencies analysis, at least the following two 
verifications must be carried out: (1) each sample or series of data 
from a site must be reviewed to look for erroneous data, that is, values 
that are too large or extremely low, as well as repeated, which could 
originate in the transcription and (2) the samples should be compared 
with each other and with the closest ones, to verify a change in 
magnitude, for example as the size of the basin grows or its location 
varies from one area to another from the analyzed region (Hosking € 
Wallis, 1997; Campos-Aranda, 2010) (Table 1). 
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Table 1. Critical values of discordance (D;) according to the number 
of stations of the group under analysis (Hosking € Wallis, 1997). 


13 2.869 


14 2.971 


Fortunately, the erroneous values, the scattered events 


(outliers), the trend and the leaps or changes in the data mean, are 
reflected in the L moments of the sample (equations (2) to (5)). 
Therefore, a convenient mixture of the L ratios in a single statistic ( D;) 
that measures the discordance between the L ratios of the site and the 
group averages has been suggested as a basic test to detect sites that 
are discordant with the group as an everything (Hosking € Wallis, 
1997). 


Assuming there are NS sites in the group being analyzed. u; = 


[tititi]” is defined as a vector containingthe sample L ratios: tz, ta, and 
ta of each site j, defined by equations (6) to (8). The superscript T 


means transposed since uj is a row vector. The average (unweighted) 
vector of the group will be (Hosking 8 Wallis, 1997): 
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The matrix A of the sum of squares and cross products will be 
defined as: 


A= 35 (4-8) (y -8) (11) 
Finally, the measure of the discordancy of each site will be: 
D)=5(y-1) 47 (4-1) (12) 


So, when D; is greaterthan the critical values (D¿) the site will be 
discordant with the group. 


Seasonality indices 


Zrinji and Burn (1996), Bum (1997), Cunderlik and Burn (2002), and 
De-Michele and Rosso (2002) have suggested that the mean date of 
flood occurrences and the seasonal regularity can be used as measures 


of similarity of the hydrological response of a basin. Therefore, such 
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basins can be integrated into a region for regional frequency analysis. 
This regionalization approach has the advantage of reserving the use 
of hydrometric information to verify regional homogeneity. In this VEA 
study, the average occurrence date of the maximum monthly runoff 
volume of each year and the seasonal regularity were evaluated. 


Directional statistics date from the early 1970s and are a simple 
tool to define measures of similarity from the dates of occurrence of 
the events analyzed. In Mexico they were introduced by Ramírez- 
Orozco, Gutiérrez-López and Ruiz-Silva (2009). 


There are various conventions or ways of using the circle to 
define directional statistics; From now on, the one that is similarto the 
clock hours will be used, since the 12 months of the year are analyzed, 


so January is one, June six, September nine, and December twelve. 


To establishthe seasonality indices, the first step is transforming 
the occurrence month of each maximum monthly runoff volume in each 
year, to an advanced angle in degrees (6;), with a 30% increase in each 
month; with 15% for January, with 75% for March, 135% for May, 2250 
for August, 285% for October, and so on. 


Taking n data from the month in which the maximum monthly 
runoff volume of each year occurred, the average coordinates of the 


mean angle (9) of the monthly maximums will be: 


x= 2)", sen(0,) (13) 


n 
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y= - 7. cos(6,) (14) 


Therefore, the mean direction 6 ofthe average occurrence of the 


months of the maximum monthly runoff in polar coordinates will be: 
E: 
0 = tan E) (15) 


The application of the previous equation is carried out by first 
obtaining the angle alpha (a) tangent of x between and y, both with a 
positive sign; then if x and y are positive 9 = a, ifx>0 and y<0 9=180"- 
a, if both are negative 9=180%+a and finally, if x<0 y y>0 9=360%-a. 

Basins with similar values of 6, can be expected to show 
similarities in otherimportant hydrological characteristics. The angle 9 


may be related to the size of the basin and its geographical location 
within the analyzed region. 


A measure of the variability of the n occurrence dates of the 
maximum monthly runoff of each year, to the 6, can be estimated by 


calculating the resulting mean, which is expressed as: 


P=/2+77 (16) 
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The seasonality index is a dimensionless measure of the data 
dispersion, ittakes values between zero and one. A unit value indicates 
that all the maximum monthly runoffs occur in the same month, while 
a value closeto zero implies great variability of occurrences throughout 
the year. 


Ordinary moments linearity test 


Gupta, Mesa, and Dawdy (1994) found that if a regionis homogeneous 
conceming the hydrological variable X;, then its ordinary moments 
(MOs) show logarithmic linearity with the basin areas (4). Therefore, 
the following expression (Ferro €: Porto, 2006) allows for verifying such 
linearity: 


MO, = 


Ea cons = 1,2,.3,4 dl 


being, s the order of the moment, n the number of elements in the 
record and X;(4) the data associated with the value of the basin area. 
Taking to a semilogarithmic graph, the basin sizes in the abscissa in 


the logarithmic scale and the ordinate, in arithmetic scale, the decimal 
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logarithms of MOs, the slopes of each grouping or cloud of moments of 
order s can be calculated, using the following type of linear regression 


lines: 


MO, =b+m:(logA) (18) 


Test of the linear regression of the B, 


Kumar, Guttarp, and Foufoula-Georgiu (1994) proposed a test that 
verifies regional hydrological homogeneity, based on an equation 
similar to the previous one, but that uses weighted probability 
moments (fr). This test was applied by De-Michele and Rosso (2002), 
and Ferro and Porto (2006). A version associated with the concept of 
Kumar et al. (1994) was used and verified by Varas (2000) in the 
U.S.A., France, England, and Chile, for the linear relations of Br, of a 
certain order and the immediate lower order. This means that the 
following linear regression is fulfilled in a hydrologically homogeneous 


region: 


B, =b+m:B,_, conr = 1,2,...,4 (19) 
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in which, b and m are the ordinates to the origin and the slope of the 
linear regression line (equations (22) and (23)), whose quality of fit is 
defined by the linear correlation coefficient (rxy): 


IB Br) (81 Br) 


» : 2 pa77 (20) 
En, _,-8,)pi-p,)] 


xy 


The regression line for a linear trend 


The dependent variable (y) is considered to be the annual runoff 
volumes VEA; in Mm? and the times or years T; is the abscissa (x), in 
this case, equal to the ¡-th value ¡. To test whether the slope (m) of 
the regression line adjusted by least-squares of the residuals is 
statistically different from zero, a test based on the Student's t 
distribution defined by the following equations is used (Ostle € 
Mensing, 1975): 
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b=VEA=m:-T (22) 


Me a 2-qp2 (23) 
t== (24) 
Sn 
$ E (25) 
m y? (71, DM? 
n 2 
Ss? E 2 ¡1 (VEA¡-VEA¡) (26) 


(n-2) 


VEA and T are the means; in the previous equation, VEA, is the 
value estimated with the regression line (Equation (21)). Sé and Sf, are 
the variances of the errors and the slope. If the calculated value t is 
greater than the critical one (tc), obtained for the student's t 
distribution with n - 2 degrees of freedom and a = 5 %, ina two-tailed 
test, the slope mis significant, that is, there isa lineartrend. The test 
has the limitation of not distinguishing between persistence and trend 
(Adeloye 8: Montaseri, 2002). To estimate the value of tc, the algorithm 
proposed by Zelen and Severo (1972) is used, with Z = 1.95996 for 
reliability (1 — a) of 95 %: 
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t, =Z+G1/v+6G2/v? + G3/v? + G4/v* (27) 
where: 


G1= (23 4+2)/4 
G2= (5Z* + 162? + 32/96) 
G3 = (327 + 197% + 177? — 15Z7)/384 


G4= (7929 + 77627 + 14827? — 19207? — 9457)/92160 


Available hydrometric information 


In Hydrological Region No. 10 (Sinaloa), Mexico there are 42 
hydrometric stations with reliable records and only 22 have their runoff 
regime unaffected by reservoirs and records older than 20 years, 
therefore, they were used. In the BANDAS system (IMTA, 2002) all the 
information on annual runoff volumes of such gauging stations was 
obtained, the general characteristics of which are presented in Table 2, 
shown in decreasing order of basin size. This information comes from 
CD-1, itis displayed on a monthly basis and is processed, as already 


indicated, in millions of cubic meters (Mm3). As expected, the records 
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of the Guamúchil and Huites hydrometric stations cover even before 
the year their respective reservoir construction (Eustaquio Buelna and 
Luis Donaldo Colosio) affected their regime. 


Table 2. General characteristics of the 22 hydrometric stations 
processed in Hydrological Region No. 10 (Sinaloa), Mexico. 


River or creek RP y (NA) VEMA 
(km2) (Mm3) 


E R. PE * 1942-1992 [ARA ad 057 | 4 ia] 818 


Santa Cruz R. San Lorenzo 1944-1987 (44) 8 919 | 1 655.043 


Guatenipa 1I *1969-2001 (33) | 8 252 | 1 563.583 


5 | Jaina R. Sinaloa *1942-1998 (57) | 8 179 | 1 372.602 


Palo Dulce R. Chinipas 1968-1982 (25) | 6 439 932.573 


Min rt 
Mii e Ele 
leal Masia Dil src rod Ds 


Naranjo A. Ocoroni * 1939-1984 (46) | 2 064 202.032 


13 | Acatitán R. Elota 1955-1994 (40) 1 884 417.784 
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R. Badiraguato | *1960-1999 (40) | 1018 | 252.900 
El Quelite R. Quelite 1961-1992 (32) 148.317 


18 | Zopilote A. Cabrera 1939-1986 (48) 83.262 
2 
2 


A. Sivajahui *1960-1983 (24 7.545 
Bamicori *1951-1983 (33) 12.672 


Description of acronyms: 


RP= years of registration processed. 
NA= number of years. 

A= basin area. 

VEMA= mean annual runoff volume. 


*Continuous records. 


In the fourth column of Table 2, the 13 records that are 
continuous have been indicated with an asterisk, that is, they do not 
present doubts in defining their span of full years. In the remaining 
nine, which are: Santa Cruz, Palo Dulce, Ixpalino, Toahayana, 
Tamazula, Acatitán, Choix, El Quelite and Zopilote, there are two or 
more missing years, generally afterthe year 1980 and therefore, there 
were difficulties in defining the period of years to process. 


Figure 1 shows the geographic location of the 22 basins that 


define the hydrometric stations of Hydrological Region No. 10 
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(Sinaloa), Mexico, citedin Table 2, which will be processed in your VEA 
registry. 
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Figure 1. Location in Hydrological Region No. 10 (Sinaloa), Mexico, 


of the basins of the 22 hydrometric stations processed. 
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Lapses of records to process 


In the definition of the recording periods of annual runoff volume to be 
processed, two work approaches can be defined: (1) as independent 
and identically distributed data records and (2) as chronological series. 
The first format complies with the requirement of frequency analysis 
and allows the selection of records interrupted by one or more years. 
The second format involves preserving the sequences of wet and dry 
years, and this is how runoff should be studied; therefore, it was 
adopted. The fourth column of Table 2 indicates the continuous 
recording lapses that will be processed. 


Once such continuous periods were defined, the missing monthly 
data for the incomplete years, as well as the missing years, were 
calculated, as detailed in the following section, to define the annual 
mean runoff volume (VEMA), which is cited in the last column of Table 
Ze 


Deduction of missing data 
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To estimate the missing monthly runoff volumes in incomplete years, 
a simple method was followed, which consisted of averaging the same 
month of the previous and subsequent year, to assign such value to 
the missing month. This is considered acceptable due to the regularity 
of the runoff season year after year. With this procedure, the records 
were completed at the stations: San Francisco, Santa Cruz, Guatenipa 
II, Jaina, Palo Dulce, La Huerta, Toahayana, Naranjo, Choix, and 


Pericos; in the one or two years. 


For the missing years, the weighted transport was applied with 
the inverse of the distance between hydrometric stations (Campos- 
Aranda, 2008), using two or more records. Such was the case in 1987 
atthe El Quelite station, basedin Ixpalino and Acatitán; as well as from 
the year 1985 of Tamazula, using Santa Cruz and Guatenipa II. 


When the weighted transport could not be applied, the missing year 
was deducted based on the mode of the annual register, estimated with 
the Mixed Gamma distribution; this was the case in 1983 in Urique IT, 


using 26 annual values. 


Regional method to be applied 
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The fourth column of Table 2 highlights the great variety of start and 
end dates of the annual runoff volume (VEA) records. The later 
highlights the impossibility of defining a common period, of adequate 
amplitude (> 30 years), of VEA records to be processed in the region 
analyzed with the 7£/ method. 


Given this situation, it was decided to apply the regional method 
of stations-years, to process the records that will form a homogeneous 


region from a hydrological approach. 


Search for regional PDF 1: Pearson type I1I (PT3) 


This probabilistic model has three parameters: Lu (location), o (scale), 
and y (shape). If y> O, the vaniation of X is: E < X <oo and we have 
the most common or widespread PT3 distribution; if y = O the Normal 
distribution occurs whose range of X is: -oo < X < 0 and when y < 0 
the reverse PT3 occurs, with an interval of X of —-o00 < X < 8 (Hosking 
8 Wallis, 1997). This PDF does not have an explicit solution, so the 
value X associated with a probability of non-exceedance (F) is 
estimated through the frequency factor (Kf), with the General equation 
of Hydrological Frequency Analysis (Chow, 1964): 
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X(F)=u+K;,:0 (28) 
the expression for Ks is (Kite, 1977): 
Ef =Z 4 (22 Dv. +2? 62D)? (22D +2 yi+ ly? (29) 


in which, y, =y/6 and Z is the standard normal deviation function of F, 


it is estimated with the following approximation (Zelen € Severo, 
1972): 


2.515517 +0.802853:w+0.010328 -w? 
gay (30) 
1+1.432788:'w+ 0.189269'w*+0.001308:'w 


being: 


w=y/In(1/F?) (31) 


when F > 0.50, F= 1 - Fis used in Equation (31) and the result of 
Equation (30) the sign is changed. 


The sign of tr, defines that of the asymmetry coefficient (y) in 
Equation (36). If 0<|z,|]< 1/3 we make z= 3r(1,)? to lead to the 
following expression (Hosking 8 Wallis, 1997): 
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2+0.1882:22+0.0442:z 
If 1/3< |r3| <1 we make z = 1—|r3| and applies: 
E DE AT (33) 
1-2.78861:z+2.56096 :22—0.77045-z 
Finally: 
4= Ay (34) 
o=2,'yr:a:T(a)/T(a+ 1/2) (35) 
y =2/ya (36) 


To evaluate the Gamma function, the Stirling approximation 


(Davis, 1972) was used: 


1 
EE EA 


AE 2022, Instituto Mexicano de 
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h = 12w 288w? 51840w?  2488320w* 


Search for regional PDF 2: Log-Normal (LN3) 


In this probabilistic model, the parameters are € (location), a (scale), 
and k (shape). In the Hosking and Wallis (1997) array of variables, the 
interval of the variable X is: —o <X<3+a/k when k > 0, -00 < X <oo 
if k=0 (Normal distribution) and ¿+a/k<X<o when k<0; 
furthermore, X is related to the random variable Z that has a standard 
Normal distribution (Equation (30)), defined as: 


1) de) whenk +0 (38) 


The shape parameter k is only a function of 7, with the following 
valid relation for |r¿| < 0.94, corresponding to |k| < 3: 


isa 2.0466534 —3.6544371:13+1.8396733:T3—0.20360244 «TÉ£ 
bg 


2 4 6 _ (39) 
1-2.01582173 :13+1.2420401 'T3-—0.21741801:75 


Tecnología y ciencias del agua, ISSN 2007-2422, 
Tecnología del Agua. O pen Access bajo la licencia CC BY- 13(4), 01-74. DOI: 10.24850/j-tyca-2022-04-01 
NC-SA 4.0 (https ://creativecommons.org/licenses/by- 

nc-sa/4.0/) 


mr Que: | y) Check for updates 


Tecnología y 


CienciaszAgua 


The other parameters are: 


E=24 -2(1-eX72) (40) 
E 
Era Eb) 


In the above equation, d(x) is the probability of non—-exceedance 
up to the value x of the Normal distribution, equal to: 


d(x)= 5, pat (42) 


where, p(x) is the ordinate of the Normal curve at x, equal to: 


=x? /2 


Pp) == (43) 


To evaluate Equation (42), the following approximation was used 
(Zelen € Severo, 1972): 


DP(x)= 1- 960 (ce, d+c,:d?+- +0: d”) (44) 
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d= —— (45) 


— 14+0.2316419:x 


c, = 0.31938153, c, = -0.356563782, C¿ = 1.781477937, c, = -1.821255978 y 
Cs = 1.330274429. 


Search for regional PDF 3: Potential Transformation 
(TP2) 


Chander, Spolia, and Kumar (1978) suggest using the Box-Cox 
transformation to normalize the data X; to z;, and then, applying 
Equation (28) with a Normal distribution and finally, obtaining the 
desired estimate X(F), with the i¡nverse transformation. The 


expressions cited are: 


z¡ = (XP? — 1)/p4 when p +0,X,>0 (46) 


z¡ =1In(X,) whenp=0,X,>0 (47) 
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Equation (46) is a general potential transformation, in which the 
logarithmic, reciprocal, and square root conversions are special cases 
since Qp is a constant ranging from -1.0 to 1.0. Obtaining q is not 
explicit, but an increase or decrease causes the same effect in the 
coefficient of asymmetry (Cs) of the z;, then its optimal value will be 
the one that takes the Cs to zero. 


Escalante-Sandoval, and Reyes-Chávez (2002) indicate that 
Equation (46) normalizes the Cs=0, but not the kurtosis coefficient (Ck 
+ 3). The correction of the kurtosis is achieved with a second 
transformation, exposed by these authors. According to Chander et al. 
(1978) in frequency analysis, the second transformation is not 
necessary. Equation (28) is now: 


Z(F)=zm>+2Z:zS (48) 


in which, zm and zs are the mean and standard deviation of the 
transformed data, and Z is estimated with Equation (30). The estimate 


sought is obtained with the expression: 


X(F) =[p (FP) + 111 (49) 


Standard error of fit 
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According to Chai and Draxler (2014), this is the most common 
indicator for contrasting probability distributions to real data, it was 
established in the mid-seventies (Kite, 1977). The expression for the 
standard error of fit (EEA, by its acronym in Spanish) is: 


1/2 


EEA = [Pltza VER VERO] 


(n—np) 


(50) 


VEA; ¡is the annual runoff volumes ordered from lowest to highest, 
whose number is n and VEA, the estimated annual runoff volumes, for 
the probability of non-exceedance (F) estimated with Equation (51) or 
Cunnane's formula, which according to Stedinger (2017) leads to 
unbiased estimates in most of the PDFs used in hydrology or 
probabilistic model that is contrasted. Finally, np is the number of fit 
parameters forthe PDF, with three forthe PT3 and LN3 and two for the 
TP2, 


_ 1-0.40 


(51) 


— n+0.20 
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Results and their analysis 


Homogeneity tests 


To make the results of a regional frequencies analysis reliable, the 
records of annual runoff volume (VEA) to process must have been 
generated by a stationary random process, which implies that it has 
not changed over time. Therefore, VEA records must consist of 


independent data, free from deterministic components. 


To verify the above, seven statistical tests were applied, one 
general, the Von Neumann Test, and six specifics: two for persistence 
(Anderson and Sneyers), two for trend (Kendall and Spearman), one 
for variability (Bartlett), and last, related to the change in the mean 
(Cramer). These tests can be consulted in WMO (1971), and Machiwal 
and Jha (2012). All the aforementioned tests were applied with a 
significance level (a) of 5 % and their results are shown in Table 3, 
indicating NH loss of homogeneity detected. 
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Table 3. Two statistical parameters and results of the seven 
homogeneity tests, inthe 22 hydrometric stations of the Hydrological 
Region No. 10 (Sinaloa), Mexico. 


Homogeneity tests 


pa 


Huites 


San Francisco 


05) 


Santa Cruz 
Guatenipa Il 
aina 


Palo Dulce 


EUA 


0 
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Zopilote 1.009| 0.031 
El Bledal 0.695| 0.039 
Pericos 1.776| 0.322 


La Tina 0.833| 0.017 


Description of acronyms: 


Cv= coefficient of variation (S/X), dimensionless. 

ríi= coef. of serial correlation of order one, dimensionless. 
CVN= Von Neumann quotient. 

TA= Anderson test. 

PS= Spearman test. 

TS= Sneyers test. 

TB= Bartlett test. 

PK= Kendall test. 


PC= Cramer test. 


It was found that the following six hydrometric stations showed 
an upward trend: Huites, Santa Cruz, Naranjo, Badiraguato, Zopilote, 
and Pericos. Huites, Naranjo and Pericos also show persistence, along 
with the following five: Guatenipa II, Jaina, La Huerta, Urique II, and 
Choix. It is observed that persistence is associated with positive values 
of the serial correlation coefficient of order one (ri) higher than 0.20, 
except in Badiraguato. 
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Discordancy test 


After defining the 22 records to be processed (Table 2), equations (2) 
to (9) were applied to obtain the ratios of L moments, shown in Table 
4, 


Table 4. L-moment ratios in the 22 hydrometric stations of the 
Hydrological Region No. 10 (Sinaloa), Mexico. 


tz ta 


2 San Francisco 0.23810 0.19358 0.12958 


Guatenipa II 0.23597 0.10004 0.11245 


5 Jaina 0.23033 0.17797 0.18571 
Palo Dulce 0.17121 0.24359 0.04599 
Ixpalino 0.21747 0.02728 0.10958 


La Huerta 0.26480 0.27554 0.15224 
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10 Urique II 0.29308 0.25481 0.15198 


13 | Acatitán 0.29624 0.20962 0.21919 
14 | Guamúchil 0.36576 0.37388 0.25532 


15 0.21305 0.18620 0.22018 


Badiraguato 0.32196 0.27868 0.20039 


17 | El Quelite 0.48463 0.59616 0.56783 
18 | Zopilote 0.49516 0,337.37 0.21933 
19 El Bledal 0.36091 0.22302 0.08793 


20 0.65688 0.64361 0.45957 


Then, based on the values in Table 4 and equations (10) to (12), 


the Discordance test was developed, which allowed to obtain the 
results shown in Table 5. The test revealed that the Palo Dulce records 
(Ds = 3.62) and El Quelite (Di7 = 4.12) exceeded the critical value of 
the Discordance of 3.00 (Table 1) and were therefore discordant. The 
Pericos record showed a high value (D2o = 2.71) and therefore, it was 


also discarded. 
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Table 5. Results of the Discordancy Test (D;) in the 22 hydrometric 
stations of the Hydrological Region No. 10 (Sinaloa), Mexico. 


Palo Dulce 


Jaina 


0.49 


*exceeds critical value. 


About the quotient zz, it is observed in Table 4, that the value of 
Palo Dulce was very low and those of El Quelite and Pericos very high, 
in contrast with their surrounding values by basin size, in the order of 
0.22 and 0.33, respectively. On the other hand, in Table 3 and relation 
to the value of Cv, the magnitude of Palo Dulce was very low and those 
of El Quelite and Pericos quite high; in other words, discordant. 


Tecnología y ciencias del agua, ISSN 2007-2422, 
13(4), 01-74. DOI: 10.24850/j-tyca-2022-04-01 


Tecnología del Agua. O pen Access bajo la licencia CC BY- 


NC-SA 4.0 (https ://creativecommons.org/licenses/by- 
nc-sa/4.0/) 


o) | y) Check for updates | 
OPEN ACCESS | 
Tecnología y 


CienciaszAgua 


Regional homogeneity 1: Seasonality indices 


For each year of the 19 records processed (Table 2), the month with 
the maximum runoff was detected. Such data were processed based 
on equations (13) to (16) and the results were obtained, which are 
summarized in Table 6 and drawn in Figure 2. A great seasonality is 
observed, since in most of the records ¡ts average direction (9) occurs 
in the second fortnight of August and the first of September, with 


angles that fluctuate between 230.90 in Bamícori (asterisk) and 2590 
in Urique II (black star). 
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Figure 2. Polar representation of the seasonality indices, inthe 19 
hydrometric stations of the Hydrological Region No. 10 (Sinaloa), 
Mexico. 


Table 6. Seasonality indices in the 19 hydrometric stations of 
Hydrological Region No. 10 (Sinaloa), Mexico. 
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(0) (7) (0) (r) 


pai 
A A 


Jaina 239.3|0.5105| 15 | Badiraguato | 251.7 | 0.7626 


Ixpalino 245.3|0.8479| 16 | Zopilote 238.8 | 0.8572 


Aaa 246.0 | 0.5544 El Bledal 232.5 | 0.8817 


Toahayana 255.0 | 0.6046 La Tina 245.1| 0.7836 


a Urique II 259.0 | 0.5474 230.9 | 0.8121 
sa ci IA A 


Regarding variability, the least dispersion is found in Santa Cruz 


(black-white box) and the greatest in Jaina (black diamond). It is also 
observed that six stations have a seasonality index of the order of 0.60 
or less, they are Huites (black box), Guatenipa II (black-white circle), 
Jaina, La Huerta (black-white rhombus), and Toahayana (black 
pentagon) and Urique II. These six stations have large basins greater 
than 4000 km2. 


Based on the results of Table 6, shown in Figure 2, itis concluded 
thatthe 19 VEA records were processed to form a homogeneous region 
from a hydrological point of view. 
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Regional homogeneity 2: Linearity of ordinary 
moments 


For each of the 19 records to be processed (Table 2), Equation (17) 
was applied to obtain their ordinary moments of order 1 to 4; the 
natural logarithm was applied to these values to obtain magnitudes 
easierto chart. Such results are shown in Table 7. 


Table 7. Natural logarithms of ordinary moments (MOs) in the 19 
hydrometric stations of Hydrological Region No. 10 (Sinaloa), Mexico. 


A [E 370 [EEE 883 pea 517 Mid 250 


IN Francisco 7. Missial 15. Asael => 908 Desa 149 


IO Ii a cd a 


Md iocecciall Mil! ol ha Fai 
Jaina 7.224 14.610 | 22.138 29.787 
MIR NO o bic Rol sia 
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Toahayana 7.005 14.138 | 21.370 28.678 


Urique II 6.075 12.400 | 18.936 25.633 


5.308 11.264 | 17.656 24.332 


12 Acatitán 6.035 12.339 | 18.884 25.642 


14 


Choix 5.686 1.13 17.475 23.567 
Badiraguato 5.533 11.376 | 17.478 23.766 
16 Zopilote 4,422 9.536 15.074 20.852 


17 El Bledal 3.821 8.029 12.562 17.378 


Ordered to the origin (b) -3.723 -6.554 | -8.837 | -10.787 


Slope (m) 2.823 5.466 8.010 10.502 


Coeff. correlation (r,y) 0.964 0.965 0.964 0.963 


In Figure 3, there is a semi-logarithmic graph, where the ordinary 


moments of Table 7 have been drawn in the ordinate axis, and a 
logarithmic scale for the basin areas is charted in the abscissa. The 
values of b and m of Equation (18) are shown in the respective column 
of s in Table 7. The internal graph of Figure 3 shows the linearity of the 
slopes (m) and the values of rxy confirm that the 19 hydrometric 


stations analyzed to form a hydrologically homogeneous region. 
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Figure 3. Relationship between ordinary moments and basin area in 
the 19 hydrometric stations of Hydrological Region No. 10 (Sinaloa), 
Mexico. 
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Regional homogeneity 3: Linear regression of the $, 


moments 


Now, for each of the 19 records processed (Table 2), Equation (9) was 
applied with r from 1 to 4, since r = O is mean and it is shown in Table 
2. Such results are shown in Table 8. 


Table 8. Weighted Probability Moments (B,) in the 19 hydrometric 
stations of Hydrological Region No. 10 (Sinaloa), Mexico. 


[EA pida 703 [ARS 151 | 1551. sal aia 210 


sra [spas 455 ema 915 O 060 isc 370 


Santa Cruz 1014.385 745.368 595.512 499.392 


Jaina 844.379 | 624.990 | 502.423 | 423.236 
IM ia 20 Hao acid rico] Ei 
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670.252 | 489.076 | 387.691 | 322.538 


AN Urique II 281.128 | 214.070 | 175.137 | 149.325 
366.742 | 265.611 | 210.191 | 174.859 
1 


[A A 149.155 120.729 | 102.637 90.002 


1 
270.775 | 205.468 | 167.983 | 143.351 
87.910 69.388 58.374 50.916 


14 Choix 178.797 131.614 105.574 88.851 


1 


[AR Badiraguato 167.162 | 128.794 | 106.354 91.385 


5 
16 Zopilote 62.245 50.686 43.298 38.132 


17 El Bledal 31.054 24.060 19.812 16.934 


Ordered to the origin 7.4878 3.2501 2.0339 1.4745 
(b) 


Slope (m) 0.6110 0.7352 0.7980 0.8362 
Coeff. Correlation (r;,y) 0.9999 0.9999 0.9999 0.9999 


The results of applying Equation (19) are shown in the respective 


column of r in Table 8, so, in the first one, there is the linear 
relationship between fB1 and fo, and inthe last column, the one between 
Ba and Baz. 


Based on the value of rxy (Equation (20)), it is observed that such 


relationships are almost perfect. Therefore, this method also guides us 
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to conclude that the 19 hydrometric stations in Table 8 form a 
homogeneous region from a hydrological perspective. 


Elimination of non—-homogeneous records 


In Table 3, the results of the statistical tests application of local or 
individual homogeneity showed that six stations or processed records 
have an upward trend, these are Huites, Santa Cruz, Naranjo, 
Badiraguato, Zopilote, and Pericos. The last one had already been 
eliminated inthe Discordancy test. Persistence was also detected in the 
first two and the following five: Guatenipa II, Jaina, La Huerta, Urique 
II, and Choix. These seven stations with persistence are analyzed in 
the next section. 


Based on equations (21) to (27), it was tested, in the six records 
with a trend, if it was significant; only at the Zopilote station was not. 
Because the trend detected was upward, it was decided to discard such 
records, since their data would tend to increase the annual runoff and 
this would increase the uncertainty associated with its estimation, 
through the regional frequency analysis. Due to the above and finally, 
15 records of annual runoff volume were then processed to define their 
regional probability distribution. 
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Regarding persistence 


In the annual runoff, persistence is a statistical property that originates 
in the occurrence of wet and dry years; Therefore, the processing 
registration periods of the stations with the highest value of the serial 
correlation coefficient of order one (ri) in Table 3 were compared: La 
Huerta (0.486), Urique II (0.381) and Guatenipall (0.336), to look for 
similarity in their lapses and in effect it occurs in the years from 1970 
to 2000. The other two stations that showed persistence were Jaina 
(0.216) and Choix (0.297), with different records both in length and in 


occurrence. 


Khaliq, Ouarda, Ondo, Gachon, and Bobée (2006) indicate that 
the validity of the flood predictions obtained through probabilistic 
analysis will be questionable if the observations of the record to be 
processed are dependent, that is, it has persistence. The foregoing is 
valid when looking for floods or predictions associated with very low 
exceedance probabilities since they are obtained at the extreme right 
of the probability distribution; This is not the case when generating 
synthetic annual runoff sequences, which cover the entire domain 


extension of the probabilistic model used (Salas, 1993; Cannarozzo et 
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al., 2009). Due to the above reasoning, it was not considered 
necessary to remove the stations with persistence. 


Selection of the regional PDF 


All the calculations described below were developed with dimensionless 
data, since the annual runoff volumes of each of the 15 records had 
been divided by their mean value or the Runoff Index (IE). The three 
PDFs recommended by Cannarozzo et al. (2009) were applied to each 
dimensionless record, which is: the Pearson type III (equations (28) to 
(37)), the Log-Normal of three fit parameters (equations (38) to (45)) 
and the potential transformation (equations (46) to (49)). The results 
of the PDF that led to the lowest standard error of fit (EEA, Equation 
(50)) have been condensed in Table 9. It is observed that the Pearson 
type III PDF did not achieve a better fit in any of the 15 records 
processed. 


Table 9. Best fits and estimated values for the indicated deciles in 
the 15 hydrometric stations of Hydrological Region No. 10 (Sinaloa), 
Mexico. 
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Probability (F) of the indicated deciles: 


0.060 | 0.528 | 0.733 | 0.919 | 1.151 | 1.591 
Guatenipa II TP2 | 0.061 | 0.504 | 0.745 | 0.946 | 1.180 | 1.576 


Jaina LN3 | 0.053 | 0.532 | 0.743 | 0.927 | 1.150 | 1.558 
E 0.035 | 0.516 | 0.789 | 0.988 | 1.194 | 1.502 


La Huerta TP2 | 0.057 | 0.514 | 0.695 | 0.873 | 1.117 | 1.663 


Toahayana TP2 |0.059 | 0.534 | 0.777 | 0.969 | 1.179 | 1.514 


Urique II dl ds 0.439 | 0.657 | 0.869 | 1.148 | 1.716 


0.034 | 0.614 | 0.808 | 0.965 | 1.143 | 1.439 
0.102 | 0.412 | 0.662 | 0.890 | 1.175 | 1.722 
0.074 | 0.371 | 0.567 | 0.781 | 1.101 | 1.903 


[A 0.062 | 0.570 | 0.761 | 0.929 | 1.135 | 1.517 


Zopilote LN3 | 0.166 | 0.115 | 0.402 | 0.713 | 1.164 | 2.199 


ido Mi id al oa a td 


Table 10 and Table 11 show in detail the results with the three 
PDFs, in the record with the highest EEA (El Bledal, worst fit) and with 
the lowest one (Tamazula, best fit). 
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Table 10. Fit parameters and values of the deciles in the El Bledal 
station (the one with the highest EEA) with the three mentioned 
PDFs. 


PT3 TP2 


Parámetro 2 o=0.677 Cs = 0.002 


Parámetro 3 y = 1.346 Ck = 2.622 


EEA (adim.) 0.185 WN 0.176 


P(X<x) TT v? J1E) 

0.20 0.438 0.445 

0.30 0.574 0.566 

0.40 0.711 0.690 
9 

[ERE A 1.020 0.983 

0.70 1.417 1.181 


0.80 1.478 1.455 


Table 11. Fit parameters and decile values in the Tamazula station 
(the one with the lowest £EA) with the three mentioned PDFs. 
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PT3 LN3 TP2 
Parámetro 1 u = 1.000 ¿ = 0.961 p = 0.295 
Parámetro 3 k = —0.246 Ck = 2.910 


EEA (adim.) 0.035 0.035 0.034 


0.20 0.719 0.721 0.722 


0.30 0.804 0.806 0.808 
0.40 0.882 0.883 0.887 
MIN 1.044 1.043 1.048 


0.70 1,139 LL 37 1.143 


Finally, Table 12 shows the results of the stations-years method 


with the three PDFs used, that is, combining the 15 dimensionless 
records, which add up to 539 data. It is concluded that the FDP Log- 
Normal leads to the lowest EEA. 


Tecnología y ciencias del agua, ISSN 2007-2422, 
Tecnología del Agua. O pen Access bajo la licencia CC BY- 13(4), 01-74. DOI: 10.24850/j-tyca-2022-04-01 
NC-SA 4.0 (https ://creativecommons.org/licenses/by- 

nc-sa/4.0/) 


o) y) Check for updates 
OPEN ACCESS e A] 


Tecnología y 


CienciaszAgua 


Table 12. Fit parameters and values of the deciles in Hydrological 
Region No. 10 (Sinaloa), Mexico, with the three mentioned PDFs, 15 
hydrometric stations, and 539 data. 


LN3 TP2 


Mi 
A 


Based on the results of Table 12, the regional growth curve 


(Equation (1)) was defined, which allowed to generate of the synthetic 
sequences of dimensionless annual runoff volume, in any site or locality 
of Hydrological Region No. 10 (Sinaloa), Mexico; this is: 
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0.497 
q(F)= 1896==7 10 — exp(0.401 :Z)] (52) 


in which, Z is the standard normal deviation as a function of the 
probability of non-exceedance F. In Figure 4, the first 500 data of the 
dimensionless annual runoff volume of 25 in 25 values have been 
drawn on log-normal probability paper, then the next 30 data from 5 
to 5, and finally, the final four values; in addition to the FDP Log- 
Normal of three fit parameters. 
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Figure 4. The fit of the Log-Normal distribution to the 539 
dimensionless annual runoff volumes of Hydrological Region No. 10 
(Sinaloa), Mexico. 
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Synthetic sequences of dimensionless VEA 


For the application of Equation (52), first, the size (ns) ofthe generated 
synthetic sequence of dimensionless annual runoff volume must be 
defined. Then, making use of an efficient random number generator 
with uniform distribution (Metcalfe, 1997), the required ns numbers u; 
are obtained; Such values are made equal to the probability of non— 
exceedance (F), and Equation (30) is applied to obtain the necessary 
Z values, which are taken to Equation (52) to obtain the magnitudes 
q(F) sought. With the above process, you can calculate any number of 


synthetic sequences of V? /1EÍ (Equation (1)) equally likely to occur. 


Scaling of synthetic sequences 


By drawing on a logarithmic paper the two final columns of Table 2, 


with the basin area (4) on the abscissa and the mean annual runoff 


Tecnología y ciencias del agua, ISSN 2007-2422, 
Tecnología del Agua. O pen Access bajo la licencia CC BY- 13(4), 01-74. DOI: 10.24850/j-tyca-2022-04-01 
NC-SA 4.0 (https ://creativecommons.org/licenses/by- 

nc-sa/4.0/) 


a Check for u dates. 
OPEN Quecess | Ln) P j 


Tecnología y 


iencias¿Agua 


volume (VEMA) on the ordinate, a straight line is defined as shown in 
Figure 5 and whose equation is: 


IE! = VEMA= 0.015324- 41281914 (53) 


In the previous expression, the correlation coefficient was 0.966, 
with the 15 pairs of data used. 
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Figure 5. Relationship between the annual mean runoff volume 
(VEMA) and the basin area (4), in the 15 homogeneous stations of 
Hydrological Region No. 10 (Sinaloa), Mexico. 
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Equation (53) is valid in the entire Hydrological Region No. 10 
(Sinaloa), Mexico, and allows estimating the value of the runoff index 
(TEN), with which dimension in Mm3 is assigned to the synthetic 


sequences of VEA generated with Equation (52). 


Gathering of main results 


From Hydrological Region No. 10 (Sinaloa), Mexico, 22 records of 
annual runoff volume (VEA) were processed. The Discordancy test 
detected two anomalous records, El Quelite and Palo Dulce; moreover, 
Pericos station was also eliminated as it showed high discordance. 


In this regional study of the VEA, the average occurrence date 
(9) of the maximum monthly runoff volume of each year and its 
regularity (7) were evaluated. In Figure 2 two groups of stations were 
defined, those of low regularity with 0.60 or less have large basins 
greater than 4000 km?2. The rest of the stations have high regularity 
between 0.78 and 0.90. However, due to their average date of 


occurrence, they form a homogeneous region. 


The other two regional homogeneity tests, the results of which 


are shown in Figure 3 and Table 8, confirm that the 19 records 
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processed belong to a homogeneous region from a statistical point of 


view. 


Before starting the search for the regional PDF, four VEA records 
were eliminated, forshowing a significant upward trend: Huites, Santa 
Cruz, Naranjo, and Badiraguato. Table 2 shows that such stations have 
different basin areas, the first two large and the other small. In 
addition, from Figure 1 it can be inferred that due to theirlocation they 
do not define a specific subregion, where it could be accepted that such 
a deterministic component exists. Therefore, they were eliminated as 
they would be influenced, by excess, the estimates of the predictions. 


The processing of the remaining 15 VEA records was carried out 
by combining such records to apply the station-year method, since a 
common data period cannot be defined with such records. The annual 
data (VEA) were divided by the mean annual runoff volume (VEMA) of 
each record, to be processed as dimensionless values. 


Three probabilistic models have tested: the Pearson type III 
distributions (PT3), the Log-Normal of three parameters (LN3), and the 
statistical technique of the potential transformation (TP2), which is 
equivalent to the Normal PDF. The LN3 distribution was adopted, as it 
led to a lower standard error of fit. 


Finally, a logarithmic equation was defined betweenthe VEMA as 
the dependent variable and the basin area of the hydrometric station 
as the independent variable. Such an equation allows the "scaling" of 


the dimensionless predictions obtained with the regional PDF. 
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Comparison with other studies 


Regional frequencies analysis has been carried out predominantly for 
floods, both in Mexico and in the world. On the other hand, the runoff 
index method (1£) described in this work has not been applied in our 
country and worldwide, the need to estimate annual runoff volumes for 
the hydrological design of reservoirs has been resolved by methods 
that simulate the rainfall-runoff process (Beven, 2001; Wagener, 
Wheater, 8 Gupta, 2004). 


The application of regional frequencies analysis to estimate the 
annual runoff volume, in Sicily, Italy, by Cannarozzo et al. (2009) 
supported the objective of this work. However, both studies are 
completely different. In Sicily, three hydrological regions were defined, 
based on the altitude and the mean annual rainfall in each basin; its 
verification was carried out through the H test for heterogeneity 
(Hosking € Wallis, 1997, Campos-Aranda, 2010). 


On the other hand, from the analysis of flood frequencies carried 
out in Sicily by Ferro and Porto (2006), the tests of Gupta et al. (1994) 
and Kumar eta/. (1994) were taken, to verify the homogeneity of the 
Hydrological Region No. 10 (Sinaloa), Mexico. 
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Conclusions 


A regional frequencies analysis was developed by processing 22 
records of annual runoff volume (VEA), from Hydrological Region No. 
10 (Sinaloa), Mexico. To select the records that form a homogeneous 
region, according to a hydrological approach, the Discordancy Test was 
applied and three records were eliminated. 


The verification of the regional hydrological homogeneity of the 
19 remaining records was carried out based on three statistical tests: 
(1) seasonality indices; (2) linearity of the ordinary moments and (3) 
linear regression of the probability-weighted moments. 


Based on seven tests of local homogeneity, another four records 
that showed an upward linear trend were eliminated. Therefore, 15 
VEA records were finally processed, with record amplitudes that ranged 


from 24 to 57 years and amounted to 539 data. 


To define the regional growth curve [q(F)], three probabilistic 
models were applied: (a) Pearson type III, (b) Log-Normal, and (c) 
Normal through the Potential Transformation. The first two were 
adjusted with the L moments method. It was found, through the 
smallest standard error of fit that the Log-Normal distribution better 


represented the 539 dimensionless annual runoff values and defined 
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the following equation, a function of Z which is the standard normal 
deviation: q(F) = 0.896 — 0.497 + (1 — exp(0.401-Z)/0.401. 


Finally, for the sca/ling ofthe syntheticsequences generated with 
the previous equation, the following empirical relationship was 
established: VEMA = 0.015324 : 4128%911 | 


For basins without hydrometry, it is essential to count with a 
method that allows estimating VEA records, of any size and equally 
likely to occur as a historical record. The additional advantage of such 
synthetic sequences, having, due to their random nature, absence of 


persistence and thus, leading to a more reliable reservoir sizing. 


Therefore, the systematic application ofthe runoff index method, 
according to the procedure outlined, is recommended in other 
hydrological regions of the country to find the corresponding regional 
equations of q(F) and VEMA = f(:). 
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